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ABSTRACT 

Explosions of type Ic supernovae (SNe Ic) are investigated using a relativistic hydrodynamic code to 
study roles of their outermost layers of the ejecta in light element nucleosynthesis through spallation 
reactions as a possible mechanism of the "primary" process. We have confirmed that the energy dis- 
tribution of the outermost layers with a mass fraction of only 0.001 % follows the empirical formula 
proposed by previous work when the explosion is furious. In such explosions, a significant fraction of 
the ejecta ( >0.1 % in mass ) have the energy greater than the threshold energy for spallation reac- 
tions. On the other hand, it is found that the outermost layers of ejecta become more energetic than 
the empirical formula would predict when the explosion energy per unit ejecta mass is smaller than 
r~> 1.3 x 10 51 ergs/ M©. As a consequence, it is necessary to numerically calculate explosions to estimate 
light element yields from SNe Ic. The usage of the empirical formula would overestimate the yields by 
a factor of £ 3 for energetic explosions such as SN 1998bw and underestimate the yields by a similar 
factor for less energetic explosions like SN 19941. The yields of light elements Li, Be, and B (LiBcB) 
from SNe Ic are estimated by solving the transfer equation of cosmic rays originated from ejecta of SNe 
Ic and compared with observations. The abundance ratios Be/O and B/O produced by each of our SNe 
Ic models are consistent with those of metal-poor stars. The total amounts of these elements estimated 
from observations indicate that energetic SNe Ic like SN 1998bw could be candidates for a production 
site of Be and B in the Galactic halo only when the fraction of this type out of all the SNe was a factor of 
> 100 higher than the value estimated from current observational data. This primary mechanism would 
predict that there are stars significantly deficient in light elements which were formed from the ISM 
not affected by SNe Ic. Since this has no support from current observations, other primary mechanisms 
such as the light element formation in superbubbles are needed for other types of SNe. The observed 
abundance pattern of all elements including heavy elements in metal-poor stars suggests that these two 
mechanisms should have supplied similar amounts of Be and B. Our calculations show that SNe Ic can 
not produce an appreciable amount of Li. 

Subject headings: nucleosynthesis — relativity — shock waves — cosmic rays — supernovae: general 



1. INTRODUCTION 

The amounts of the light elements Li, Be, and B (LiBcB) 
at present are thought to be the sum of the products of the 
big bang nucleosynthesis and the subsequent cosmic-ray 
spallation reactions. The contribution from the big bang 
nucleosynthesis to 7 Li is thought to be most pronounced 
among these three elements. Though some fractions of 
Li have been depleted inside some cool stars, its abun- 
dances on the surfaces of metal-poor stars are the only 
information that we can use to identify the contribution 
from the big bang nucleosynthesis (e.g., Ryan, Norris, & 
Beers 1999). Thus it is important to know the contribution 
from the cosmic-ray spallation reactions to Li in the early 
stages of the Galaxy evolution during which metal-poor 
stars were formed. With this knowledge, we can constrain 
the cosmological parameters to synthesize light elements 
in the big bang nucleosynthesis (Ryan et al. 2000). 

In addition, the investigation of the evolution of these 
elements in the early stages of the Galaxy enables us to un- 
derstand the origin of cosmic-rays responsible for the light 
element nucleosynthesis. Recent observations of extremely 
metal-poor stars apparently suggest that the primary (not 
secondary) spallation process dominated the light element 
nucleosynthesis in the early Galaxy (Duncan, Lambert, 



& Lemke 1992). In other words, cosmic-rays composed 
of C/O interacting with protons and/or He nuclei in the 
interstellar medium (ISM) have predominantly produced 
LiBeB. A primary mechanism, in which LiBeB were pro- 
duced in supernova ejecta-enriched superbubbles, was sug- 
gested to explain the Be evolution in the early Galaxy 
(Higdon, Lingenfelter, & Ramaty 1998). Since then, two- 
component models in which light elements are produced 
by standard Galactic cosmic-rays and metal-enriched par- 
ticles in superbubbles have been investigated by several 
authors (Ramaty et al. 2000; Fields, Olive, Vangioni-Flam, 
& Casse 2000; Suzuki & Yoshii 2001). These studies con- 
cluded that a primary mechanism is needed to explain the 
observed abundance trends of Be and B with O and Fe. 
Recently, Suzuki & Yoshii (2001) investigated the chemical 
evolution of LiBcB in the early epoch of the Galaxy us- 
ing an inhomogeneous chemical evolution model developed 
by Tsujimoto, Shigeyama, & Yoshii (1999) in which star 
formation is assumed to be induced by supernova (SN) ex- 
plosions. Their model succeeded in reproducing not only 
the observed metallicity distribution of Galactic halo stars 
but also the observed abundance correlations of heavy el- 
ements. Suzuki & Yoshii (2001) also considered two ori- 
gins of Galactic cosmic-rays that synthesize light elements. 
One is from the ISM accelerated by SN shocks, and the 
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other from SN ejecta accelerated by the SN shock. The 
authors claimed that ^2 % of Galactic cosmic-rays must 
be originated from the SN ejecta to reproduce the behavior 
of the abundance of Be at the metal-poor ends. However, 
their model suffers from shortage of energy supply from 
each SN to cosmic-rays. Moreover, the energy distribu- 
tions of the cosmic-rays from these two origins might be 
different, though Suzuki & Yoshii (2001) assumed they are 
the same. 

Fields (1996) found from a numerical model of type Ic 
supernovae (SNe Ic) (Nomoto, Filippenko, & Shigeyama 
1990) that the outermost C/O layers of an SN could at- 
tain energies beyond the threshold value to produce LiBcB 
(~30 MeV per nucleon for O+H^Be, which corresponds 
to the Lorentz factor of ^1.03). Energetic explosions of 
massive stars with stripped H-rich and He layers might 
be able to produce cosmic-rays enriched with C/O in the 
outermost layers of ejecta. In the model of Nomoto, Fil- 
ippenko, & Shigeyama (1990), the outermost layers of the 
ejecta do not become so relativistic. It may, however, be 
attributed to their coarse zoning in the outermost layers, 
which has to be as accurate as possible for this purpose. 
Furthermore, their numerical calculations were performed 
with a hydrodynamic code that does not take into account 
relativistic effects. Thus it was impossible to accurately es- 
timate the contribution from SNe Ic to the light element 
nucleosynthesis with their results. 

The phenomena taking place when the supernova blast 
wave hits the surface of a relatively compact star were 
investigated in a more sophisticated fashion by Ensman 
& Burrows (1992) and Blinnikov et al. (2000) with their 
radiation-hydrodynamic codes. Their codes allow the ra- 
diation and the gas to go out of equilibrium. Both of 
their codes were non-relativistic except that they included 
light-travel-time corrections. They were concerned with 
the shock breakout of SN 1987A, because the most de- 
tailed observations immediately after the shock breakout 
were available for this SN. One of their main objectives 
was a detailed modelling of the UV burst immediately 
after the shock emergence. In addition to these numer- 
ical approaches, there have been semi-analytic approaches 
to the shock emergence in the plane-parallel medium in 
which the flow is assumed to be self-similar (Gandel'man 
& Frank-Kamenetskii 1956; Sakurai 1960). Kazhdan & 
Murzina (1992) took into account the sphericity of the 
flow in a neighborhood of the surface. These semi-analytic 
approaches are also limited to the non-relativistic flow. 

After the emergence of a very bright type Ic super- 
nova, SN 1998bw, was found to be associated with a 7-ray 
burst GRB 980425 (Galama et al 1998) and Kulkarni et al. 
(1998) inferred from their radio observations that the radio 
shock front of SN 1998bw was moving at relativistic speeds 
with the Lorentz factor between 1.6 and 2, the relativistic 
motion from supernova explosions has been investigated. 
Matzner & McKcc (1999) estimated how much mass of 
relativistic ejecta could be obtained from the explosion 
of a massive star and derive an empirical formula giving 
the mass of relativistic ejecta from explosions of stars with 
simple polytropic density structures. Later, Tan, Matzner, 
& McKee (2001) revised the empirical formula using their 
relativistic hydrodynamic code. Their formula would give 
~ 10~ 6 Mq of relativistic ejecta (Lorentz factor > 2) from 



the explosion of a star with a mass of 10 Mq and an energy 
of 10 52 ergs. 

Using the empirical formula of Matzner & McKee (1999) 
for the energy distribution of SN ejecta, Fields et al. (2002) 
concluded that SNe Ic, especially energetic events like SN 
1998bw significantly contribute to light element nucleosyn- 
thesis through spallation reactions. The energy distribu- 
tion of particles in SN ejecta used in Fields et al. (2002) 
follow a power law with a power index of —4.6. This power 
index is quite different from that of the energy distribu- 
tion of the ISM accelerated by SN remnant shock (e.g. 
Meneguzzi, Audouze, & Reeves 1971). Then Fields et 
al. (2002) adopted a "thick target" approximation instead 
of solving the cosmic-ray transfer equation to derive the 
yields of LiBcB from the energy distribution of SN ejecta. 

There still remain some problems to be addressed in 
their work; the model for stellar structures Matzner & 
McKee (1999) assumed is so simple that it involves sus- 
picions of inaccurate estimates. In addition, they set the 
adiabatic index 7=constant over the whole stages of explo- 
sions, which may vary according to the physical conditions 
and affect the resultant energy distribution of the ejecta. 

Our objective is to construct a realistic model for SN 
ejecta moving at relativistic speeds and to investigate their 
contribution to light element syntheses. To improve the 
above situations, we calculate the energy distributions of 
ejecta as a result of SN explosions using a relativistic hy- 
drodynamic (RHD) code with realistic numerical models 
for massive stars as the initial conditions. The atmo- 
spheres in radiative equilibrium are attached to the orig- 
inal models for massive stars (see references in Table 3). 
This is essential to investigate the energy distribution of 
ejecta at highest energies after the shock breakout. We 
also verify the validity of 7=constant by using more re- 
alistic equation of states that incorporates radiation and 
ideal gas in thermodynamic equilibrium. To investigate 
the change of the energy distribution of relativistic ejecta 
transferring in the ISM, we solve the transfer equation that 
takes into account the energy loss due to ionization and 
spallation reactions. From these calculations, the amounts 
of synthesized LiBcB through primary spallation reactions 
are obtained. 

In the next two sections, we describe our RHD code 
in Lagrangian coordinates (§2) and initial conditions (§3). 
Then, we show the results in §4, compare them with that 
of other authors, and discuss the differences. In §5, we 
estimate the yields of LiBeB using our explosion models 
together with the leaky box model. 

2. NUMERICAL CODE FOR RELATIVISTIC 
HYDRODYNAMICS 

We have constructed a numerical code for RHD in La- 
grangian coordinates based on the formalism presented in 
Marti & Miiller (1996). Our RHD code well reproduces 
the exact solutions of test problems presented in Marti & 
Miiller (1996), such as relativistic shock tubes. 

2.1. Equations 

In Eulerian coordinates, equations of RHD of a perfect 
fluid can be written in the following form: 

^+V.(£v) = 0, (1) 
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f + V(S.v + p) = 
| + V.(S-Dv), 



(2) 



(3) 



where t, D, S, r denote the time, the rest-mass density, the 
momentum density, and the energy density in a fixed frame 
where the fluid moves at speed v, respectively. These vari- 
ables are related to quantities in the local rest frame of the 
fluid through 

D = P W (4) 



S = phW\ 
t = phW 2 - p-D, 



(5) 
(6) 



where p, p, W and h denote the proper rest-mass den- 
sity, the pressure, the fluid Lorentz factor, and the specific 
enthalpy, respectively. The specific enthalpy is given by 



h=l + e+-, 
P 



(7) 



where s is the specific internal energy. 

We can describe these equations with a Lagrangian 
coordinate m. Since we are concerned with spherically 
symmetric explosions, we will rewrite them in the 1- 
dimensional spherical coordinate system. 



dV i7r d{r 2 v) 
dt dm 



0. 



h 47rr" 

at 



dp + GM r 
dm r 2 



dQ | ^ d(r 2 pv) = Q 
dt dm 



(8) 
(9) 

(10) 



where v denotes the radial velocity, G the gravitational 
constant, M r the mass included within the radius r, V = 
1/.D, s = VS ■ r/r (r: the radial vector) and 

0-rV-°* ' (11) 
r 

The Lagrangian coordinate m is related to r through 

47rr 2 

dm = —dr. (12) 

The gravity is included in the weak limit where no general 
relativistic effect is prominent. The pressure is described 
as the sum of the radiation pressure and the gas pressure: 

kpT 



pm p 



(13) 



where a denotes the radiation constant, T the tempera- 
ture, k the Boltzmann constant, p the mean molecular 
weight, and m p the mass of proton. The specific internal 
energy e is given by 

aT A 3kT 

s= + - . (14) 

p 2pm p 

Here we will introduce two kinds of adiabatic indices, 71 
and 72 defined by 

' dlnp\ 



7i 



72 



d In p ) s ' 
V 



1. 



pe 



(15) 
(16) 



The difference scheme of equations (8)-(10) can be writ- 
ten in the following form; 





-(r 2 v)^ 1/2 + (r 2 v) i+1/2 
(r 2 p)i-i/2 - (r 2 p)i+i/ 2 
n2 pv)i-i/ 2 - {r 2 pv) l+ i/ 2j 



(17) 

We use the Godunov method to numerically solve equa- 
tions (17) distinguishing these two adiabatic indices in our 
Riemann solver. 

The number of zones are listed in Table 3. To obtain ac- 
curate energy distributions of SN ejecta above the thresh- 
old energies for cosmic-ray spallation reactions, our models 
have zones in the outermost layers with masses well below 

icr 9 Mq. 

2.2. Boundary Conditions 

When wc calculate the explosion of a star, we need to 
set the boundary conditions at the outer edge of the star 
as follows; 

Pimax+1/2 = 0, 

Vimax+l/2 — Vimaxj 



(18) 



where imax is the zone number corresponding to the out- 
ermost layer. 

At the center, which corresponds to the zone interface 
i = 1/2, we consider an imaginary zone i — 0, where the 
physical quantities are given by 



Pa =Pi, Pa 
v = -vi, r - 



Pi 

-n, 



(19) 



to keep the symmetry with respect to the center. 

3. INITIAL CONDITIONS 

We consider four stellar models immediately before the 
core collapse as the initial conditions. These stars are orig- 
inated from 12— ~ 40 Mq main-sequence stars. Three of 
them are thought to have undergone intense stellar winds 
and lost their H-rich and He envelopes. As a result, the 
stellar surfaces mainly consist of carbon and oxygen at ex- 
plosion. These three models have corresponding real type 
Ic supernovae as shown in Table 3. These stars have be- 
come fairly compact with radii less than the solar radius. 
This compactness results in higher pressures at the shock 
breakout compared with explosion of a star with an ex- 
tended envelope if the explosion energies per unit ejecta 
mass are the same. The other star is thought to be the pro- 
genitor of SN 1987A, one of the best studied supernovae. 
We have calculated the explosion for this supernova and 
compare the result with the previous work to check our 
numerical code. 

To initiate explosions, we first replace the central Fe core 
with the point mass located at the center, release the en- 
ergy in the innermost several zones in the form of thermal 
energy (or pressure), and trace the evolution of physical 
quantities. The calculations are stopped when the ejecta 
keep expanding homologously. The parameters of each 
model are tabulated in Table 3. 
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Table 1 

Parameters of stellar models. 



M ms a 




i?* c 




M G d 




Number of zones 


SN C 


Reference 


(M Q ) 


(M ) 


(Re) 


(M e ) 


(M ) 


(xl0 51 ergs) 








-40 


15 


0.29 


13 


10 


30 


389 


1998bw 


Nakamura et al. (2001) 


20-25 


4.6 


0.32 


2.9 


2.1 


4 


398 


2002ap 


Mazzali et al. (2002) 


12-15 


2.4 


0.24 


0.99 


0.66 


1 


294 


19941 


Iwamoto et al. (1994) 



"The stellar mass on the main sequence. 
6 The stellar mass at core collapse 
c The stellar radius at core collapse 
d The mass of oxygen in the ejecta 
e Corresponding real supernova 



4. THE ENERGY DISTRIBUTION OF SUPERNOVA EJECTA 

4.1. SN 1987A 

To check our numerical code, we have performed a cal- 
culation for a type II SN 1987A and compare the resultant 
velocity distribution as a function of enclosed mass with 
that calculated by Shigeyama & Nomoto (1990) with a 
non-relativistic hydrodynamic code. Results of these two 
calculations are found to be in agreement to a fairly good 
accuracy ( ^ 0.2 % on average). 

4.2. SN 1998bw 

The explosion energy has been found to be considerably 
higher than 10 51 ergs from modelling the observed light 
curve and spectra of SN 1998bw. Asphericity has been 
suggested from the observed feature of emission lines of O 
and Fe (Maeda et al. 2002). The fact that spherical models 
cannot explain the light curve both in the early and late 
phases also suggests an aspherical explosion : The spheri- 
cal explosion with an energy of ~ 5 x 10 52 ergs can explain 
the early phases but the light curve was predicted to de- 
cline too fast to be compatible with the observations, while 
that with a lower energy ~7x 10 51 ergs cannot evolve as 
fast as observed in the early phase but can explain the 
observed decline of the late light curve (Nakamura et al. 
2001). 

Here we have calculated an explosion of a <~ 15 Mq 
C/O star with an explosion energy E ex of 3 x 10 52 ergs for 
this energetic supernova (Table 3). The mass of the ejecta 
M e j becomes 13 Mq with the rest of ~ 2 M Q collapsed 
to become a compact object at the center. The resultant 
energy distributions of the ejecta are plotted in Figure 1 
with that predicted from the empirical formula proposed 
by Fields et al. (2002): 

M>A =2 2 x 1Q -4 ( E cx \ 3 6 / M^_\ - 3 6 
M cj ' v 10" ergs J \l M e J (2Q) 

x ( 6 V 3 6 
X UOMeW 

where M(> e) denotes the mass of ejecta that have parti- 
cle energy per nucleon greater than e. At the high energy 
tail, the energy distribution from our calculation with re- 
alistic equation of states (the dashed line in Fig. 1) show 



the same dependence on energies as the above empirical 
formula but with somewhat (~ 13 %) smaller amounts of 
matter. The adiabatic indices (71 and 72) fixed to 4/3 
also agree with the empirical formula (the dotted line in 
Fig. 1). Thus this agreement is ascribed to the features of 
this particular model that the density distribution of the 
star is well approximated by a simple formula presented 
in Matzner & McKee (1999) and that the adiabatic in- 
dices in the shocked ejecta are close to 4/3. At around 
the threshold energies for spallation reactions (~ 5-~ 30 
MeV/nuclcon), the energy distribution from our calcula- 
tion significantly deviates from equation (20) (see Fig. 1). 
Therefore, it is doubtful to use equation (20) to estimate 
yields of light elements synthesized by the spallation reac- 
tions. The explosion of this supernova is so furious that 
the assumption of M(> e)/M e j << 1 already breaks down 
in this energy region. In summary, equation (20) always 
overestimates the mass with given e. In particular, it over- 
estimates the mass with e greater than the threshold en- 
ergies by a factor of ;C 3. 

4.3. Other type Ic supernovae 

The other type Ic supernovae listed in Table 3 are explo- 
sions with explosion energies per unit mass smaller than 
SN 1998bw model. The energy distribution of the ejecta of 
these supernovae follow a power law with the same index 
in equation (20) but with a significantly greater amount 
of matter accelerated to given energies per nucleon. A low 
E ex /M c j leads to adiabatic indices in the shocked matter 
greater than 4/3. The stiffer equation of states leads to a 
more efficient acceleration. The same is true for a model 
of SN 1998bw with a less explosion energy (e.g. 3 x 10 51 
ergs) as shown in Figures 2 and 3. As a result, equation 
(20) underestimates the mass of ejecta with given e for 
these supernovae. 

4.4. Energy distribution with a variable adiabatic index 

The energy distributions (— dM(> e)/de) for all the 
models have smaller power-law indices (~ —4.6) than that 
of the cosmic-rays accelerated through the shock-fronts of 
SN remnants ~ —2.6 (e.g. Meneguzzi, Audouze, & Reeves 
1971). This will affect the ratios of light elements synthe- 
sized from the cosmic-ray spallation reactions (see §5.2). 
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Fig. 1. — The integrated energy distribution M(> e) of ejecta for SN 1998bw model with _E ox = 3 X 10 52 ergs. The solid line represents the 
fitting formula of Fields et al. (2002). The results of our numerical calculation with the constant adiabatic index (71 = 72 =4/3) is shown by 
the dotted line. The dashed line shows the result with the variable adiabatic index. 




From the discussion in the preceding sections, it follows 
that equation (20) holds when the adiabatic index is very 
close to 4/3 and M(> e)/M e j ^ 10~ 5 . In particular, our 
calculations show that the energy distribution docs not 
become a simple power-law for M(> e)/M e j £ 10~ 5 . This 
must be due to the fact that the simple stellar density dis- 
tribution Matzner & McKee (1999) assumed can reproduce 
the density distributions in realistic stellar model only for 



the outer 0.001 % layers in mass. When the explosion 
energy is furious, the adiabatic index is close to 4/3 but 
the mass with e greater than the threshold energies for 
spallation reactions can become significantly greater than 
10~ 5 M cj as in SN 1998bw model. Thus equation (20) does 
not give a good estimate for the mass of ejecta in the re- 
gion between e and e + de around the threshold energies 
(-dM(> e)/de) as well as M(> e). 
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Fig. 3. — The distributions of adiabatic index 71 as a function of the enclosed mass (including the central remnant mass) immediately after 
the shock breakout for the SN 1998bw models. The solid line shows the model with _B CX = 3 X 10 51 ergs. The dashed line with E cx = 3 X 10 52 
ergs. 



The smaller the explosion energy per unit mass E ex /M e j 
becomes, the more significant our results deviate from 
equation (20) even in the high energy tails (Fig. 2) . Our 
calculations for all the models listed in Table 3 with a 
constant adiabatic index equal to 4/3 always result in the 
energy distributions of ejecta in good agreement with their 
fitting formula at high energy tails irrespective of values 
of E ex /M c j. Therefore this deviation should be caused by 
variable adiabatic indices. At the same time, the agree- 
ment of these results with equation (20) suggests that the 
simple density distributions used in Matzner & McKee 
(1999) are good approximations of realistic stars in the 
outermost layers (0.001 % of the ejecta mass). 

From our calculations with a realistic equation of state 
(Eq. (13)) varying explosion energies for three SN models 
in Table 3, we have plotted M(> e) x e 3S /M c - 3 at e = 100 
Mev/nucleon as a function of E ex /M e j in Figure 4. For 
less energetic explosions E ex /M c j £ 6x 10 50 ergs/ M , the 
lines of the three models deviate from each other. Thus a 
single fitting formula such as equation (20) cannot account 
for all the realistic explosions. We need three different pa- 
rameters A for these three SN models to express the energy 
distribution of ejecta as 



M(> e) 



10 51 ergs 
(lOMev) 



3.4 



-3.6 



1 M G 



-3.4 



(21) 



. 10 MeV 

that takes into account variable adiabatic indices origi- 
nated from the equation of states (Eq. (13)). The constant 
A is equal to 1.9 x 10~ 4 for SN 1998bw model, 2.1 x 10~ 4 
for SN 2002ap, and 2.8 x 10~ 4 for SN 19941. 

5. LIGHT ELEMENT NUCLEOSYNTHESIS 

To investigate the role of supernova ejecta in light ele- 
ment nucleosynthesis, the modification of the energy dis- 



tribution of ejecta when they transfer in the ISM need 
to be considered because the cross sections of spallation 
reactions are sensitive to energies of particles. 

5.1. Transfer equation 

Energetic ejecta accelerated by a supernova explosion 
lose energy when they collide with neutral atoms in the 
ISM and ionize them. The ionization energy loss rate of 
element i in the ejecta through H gas, (MeV/s), is given 
by the formula (Schlickeiser 2002): 



Wi =Z 2 cil t x 1.82 x 10" 13 nHiMi 

x {1 + 0.0185 \nf3H{(3- 0.01)} 



2/3 2 



(22) 



10- 



+ 2/?3 ' 

where tiri is the number density of neutral H in the ISM, 
[3 = v/c and H(x) denotes the Heaviside step function. 
The effective charge Z c ff,i is expressed as 

Z eS ,i = Zi{l - 1.034exp(-137/3Z-°- 688 )}, (23) 

where Zi is the atomic number of the element i in the 
ejecta. 

We use the leaky-box model (Meneguzzi, Audouze, & 
Reeves 1971) and the transfer equation for the mass of the 
clement i with an energy per nucleon e at time t, Fi(e,t), 
is expressed as 

8Fi(e,t) 8[wi(e)Fi(e,t)] 



dt 



de 

Fi(e,t) , Fi(e,t) 



A e 



PVi{e), 



(24) 



where A's are the loss lengths in g cm~ 2 , p denotes the 
mass density of the ISM, fj(e) the velocity of the element 
i with an energy per nucleon of e. A osc denotes the range 
before escaping from a given system (we assume A esc = 100 
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g cm^ 2 following Suzuki & Yoshii (2001)), and A n> j 
to spallation reactions. The latter is expressed as 

UpTTlp + Tl Hc mHc 



due 



A„,i(e) = 



(25) 



Here the total cross sections of spallation reactions be- 
tween particle i and proton or He are denoted by cr Pii or 
<7H e ,i, respectively. The mass of He is denoted by mn c - 
The number densities of proton and He in the ISM have 
been introduced as n p and nn e and we assume that the 
ISM is uniform and neutral: n p = rim = 1 cm~ 3 and 
"■He = 0.1 cm~ 3 . The ionization energy loss rate is pro- 
portional to rim- The energy loss rate is, however, not so 
sensitive to the ionization state of H, since the energy loss 
rate through Coulomb collisions with free electrons in the 
ionized medium is comparable to that in a neutral medium 
and have a similar energy dependence. 

Equations (24) for C and O are numerically solved with 
initial conditions given by the energy distribution pre- 
sented in the previous section: 

F i (e,0) = -X i (e) d -^±. (26) 

Here Xi (e) denotes the mass fraction of the element i with 
an energy per nucleon of e. Equations (24) are implicitly 
integrated with respect to time. 

Through the ionization energy loss, the energy distribu- 
tion of ejecta becomes harder as time passes as shown in 
Figure 5. The effect of A osc is negligible in this calculation. 
Ionization reduce the energies of most particles below the 
threshold energies for the spallation reactions before the 
particles escape from the system. 

Since Fields et al. (2002) used the "thick target" approx- 
imation to calculate yields of light elements synthesized 
by cosmic-ray spallation reactions, we compare the energy 
distribution calculated from equation (24) with that ob- 
tained from the "thick target" approximation (Ramaty & 



Lingcnfcltcr 1975). This approximation assumes the time- 
integrated energy distribution of cosmic-rays as 



Fi(e,t)dt = 



R l (e)X l (e) M(> e) 
pVi(e) 



(27) 



where the range of element i, Ri(e), is defined as ii,(e) = 
pvi(e)e/u>i. This relation can be derived from the integra- 
tion of equation (24) with respect to time t and energy 
e omitting the terms including A's. The right hand side 
of this equation is plotted in Figure 6 with the dashed 
lines for C and O together with the time-integrated energy 
distribution obtained by solving equation (24) (the solid 
lines). It is clear that this approximation gives a fairly 
good estimate for the time-integrated energy distribution 
at around threshold energies (several to 30 MeV per nu- 
cleon). However, Fields et al. (2002) adopted a somewhat 
different approximation in which 



Jo 



Fi(e,t)dt = - 



Ri(e)Xi(e) dM(> e) 



(28) 



This will lead to a factor of ' 
of light elements. 



5.2. 



pvi(e) de 
3.6 greater amount of yields 



Yields of light elements 

Light elements (LiBcB) produced through the spallation 
reactions; C, O + H, He -^ 6,7 Li, 9 Be, 10 ' n B are considered 
here. The yield of a light clement I via the i + j — > I -\ — ■ 
reaction is estimated by the following formula; 



dm 

dt 



=nj L 



<j^^r<^ ( 29 ) 



where Ni is the number of the produced light element I, 
rij is the number density of interstellar element j, a\ j the 
cross section of i + j — > I + ■ • ■ reaction, Ai the mass 




Eex/Mej 

Fig. 4. — M(> to) X eQ- 6 /M e j as a function of E e ^/M e j where eo = 100 MeV/nucleon. M(> eo) and M e j are in units of Mq, and E ex in 
10 51 ergs. The solid, dashed, and dotted lines show models for SN 1998bw, SN 2002ap, and SN 19941, respectively. The values of M(> e) X e 3 6 
for e > 10 MeV/nucleon are almost constant for different models with the same value of E ex /M e j as long as E ex /M e j > 6 X 10 50 ergs/ Mq. 
The filled circles correspond to E ex shown in Table 3. 
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1 10 100 1000 10 4 

Energy [MeV/nucleon] 

Fig. 5. — The attenuation of the energy distribution of O for SN 1998bw model with E cx = 3 X 10 52 ergs is shown. The solid line is the 
initial energy distribution of the ejecta in the homologous expansion phase. The dashed line shows the energy distribution attenuated by 
ionizations at t = 1 Myr and the dotted line at t = 5 Myr. 



number of the element i. The integration with respect to 
energy is carried out between the minimum and maximum 
energies per nucleon in the ejecta obtained from our nu- 
merical calculations, which are denoted as e m i n and e max , 
respectively. We use the empirical formula given by Read 
& Viola (1984) for the cross sections. 

6. RESULTS AND COMPARISON WITH OBSERVATIONS 

6.1. Abundance ratios 

Table 6.2 summarizes the results. Since the ratios of el- 
ements are determined solely by the energy distribution of 



the injected materials and the energy distributions of all 
the calculated models have a similar shape, the models give 
similar isotopic ratios: 7 Li/ 6 Li~1.2 and n B/ 10 B- 2.8. 
The B isotopic ratio observed in meteorites is 4.05±0.05 
(Shima 1962; Chaussidon & Robert 1995). Thus neutrino 
spallation processes in core-collapse SNe may contribute 
to the enhancement of 11 B . On the other hand, this ratio 
observed in the local ISM is around 3.4 ± 0.7 (Lambert 
et al. 1998), which can be marginally reproduced by the 
mechanism discussed here. 

The element ratio B/Bc^ 15 from our models can ex- 




1000 



Energy [MeV/nucleon] 



Fig. 6. — The solid lines show the time integrated energy distributions of C and O in the ejecta derived from direct integration of the 
transfer equation (24). The dashed lines show the same quantities obtained with the "thick target" approximation (Eq. (27)). 
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plain this ratio on the surfaces of some halo stars because 
the average ratio of halo stars is thought to be 17—21 
and some halo stars like HD94028 show fairly low ratio as 
~ 12. The solar value 29.6 suggests that at least B in the 
solar system must have been synthesized in different sites 
in addition to neutrino processes. 

6.2. The amounts of Be and B in the Galactic halo 

Since the big bang nucleosynthesis is not expected to 
produce observable amounts of B or Be, it will be easier 
to identify which primary process is responsible for the 
enrichment of these elements as compared with Li. The 
element abundances of Be are determined for more stars 
than the other light elements Li and B. The observed mass 
ratios of Ab c /Ao on the surfaces of metal-poor stars range 
from - 1CT 9 to - 2 x 1(T 8 (Boesgaard et al. 1999). The 
result presented in Table 6.2 is in good agreement with 
this observation. This is in contrast with the result of 
Fields et al. (2002). They obtained X Bc /X ~ 10~ 6 for 
a furious SN Ic like SN 1998bw. This large amount of 
Be must be diluted by a large number of other SNe that 
do not supply such a large amount of Be to be reconciled 
with the observed Xbc/Xo ratios. They found that the 
fraction of furious SNe Ic needs to be less than 8 x 10~ 4 
out of all the SNe. Ramaty et al. (2000) had proposed 
another primary process in which SNe in superbubblcs 
yield X-q b /Xo ~3x 10~ 8 on average. As a consequence, 
the overproduction of Be from furious SNe ironically made 
Fields et al. (2002) to suggest that the supply of Be from 
furious SNe needs to be a tiny fraction of that from SNe 
in superbubbles. 

On the other hand, our result does not need to dilute Be 
with other SNe. Instead, the Be/O ratio on the surface of 
a metal-poor star need to be determined by a single SN as 
was suggested for heavy elements (Shigeyama & Tsujimoto 
1998). This might be possible if the tangled magnetic fields 
in the ISM trap accelerated SN ejecta inside the region 
where the SNR shock eventually swept up. This mecha- 
nism predicts that there must be a considerable number of 
metal-poor stars with essentially no light elements. These 
are the descendants of SNe other than SNe Ic. It is, how- 
ever, likely that this conjecture has been ruled out by cur- 
rent observations. No star without light elements has been 
reported. There is also a defect in the argument of Fields 
et al. (2002). The abundance pattern of heavy elements de- 
termined by multiple SNe in superbubbles must be at odds 
with the observations that show a large scatter in the stel- 
lar abundance ratios especially for Fe-group elements and 
r-process elements (e.g., McWilliam et al. 1995; Cayrel et 
al. 2004). To reproduce the trends of all the elements ob- 
served for metal-poor stars, both of these two sites should 
have supplied comparable amounts of Be. The quantita- 
tive discussion needs a detailed chemical evolution model 
not only for light elements but for heavy elements, which 
is beyond the scope of this paper. 

The observed mass ratio X-Qe/X-a has a value ~ 10~ n 
at [Fe/H]~ —1.5 (Boesgaard et al. 1999) which is the max- 
imum value of [Fe/H] attained in the Galactic halo. The 
Fe yield from core-collapse SNe of 2.1 x Xp e & per unit SN 
mass (Tsujimoto et al. 1997) indicates that the yield of Be 
per unit supernova mass is 



in the halo. The corresponding values for the three models 
in Table 6.2 are - 2.5 x 10~ 9 (SN 1998bw), - 1.7 x 10~ 10 
(SN 2002ap), and - 4.8 x 10~ n (SN 19941). Observations 
suggest that the fraction of furious SNe Ic is ~ 1/700 of 
core-collapse SNe (Podsiadlowski et al. 2004). This frac- 
tion should have been a factor of > 100 higher in the early 
stage of the Galaxy for furious SNe Ic to have been a ma- 
jor production site for Be. The element abundance ratios 
B/Bc from the present models consistent with those for 
Galactic halo stars indicate that the same is true for B in 
the halo. 

6.3. Li in the Galactic halo 

The observed number ratios Li/H on the surfaces of 
metal-poor stars suggest that the big bang nucleosynthe- 
sis leads to Li/H^ 1.23^o 32 x 10~ 10 and have a significant 
correlation with the metallicity (Ryan et al. 2000). If this 
correlation is extrapolated to [Fe/H] = —1.5, then the ratio 
Li/H becomes ~ 4.6 x 10 -10 . Therefore the amount of Li 
produced in the Galactic halo was ~ 2.5 x 10~ 9 Mh, where 
Mh denotes the mass of H. Thus an argument similar to 
the preceding section results in the yield of Li per unit 
supernova mass of 

MLi : (31) 



Mi 



1.2 x 10" 



SN 



M SN 



5 x lO" 10 , 



(30) 



The corresponding value of SN 1998bw model is 1.4 x 10~ 8 
and the other models give less values. If this value in equa- 
tion (31) is taken at face value, SNe Ic cannot produce this 
amount of Li. However, there seem to be few stars with 
the metallicity range of —2 ^ [Fc/H] — 1 for which the 
Li abundance is measured (Ryan et al. 2000). We need to 
wait for abundance information in this metallicity range 
to deduce a firm conclusion. 

At the low metallicity end ([Fe/H] ~ —4), the observed 
Li/O ratios ~ 10~ 4 (Ryan et al. 2000) are more than two 
orders of magnitude greater than that in our SN 1998bw 
model (see Table 6.2). The Li produced from this pri- 
mary mechanism is likely to be negligible compared to the 
primordial Li. 

7. CONCLUSIONS AND DISCUSSION 

We have performed numerical calculations for SNe Ic 
explosions using a relativistic hydrodynamic code to in- 
vestigate how much mass of ejecta is accelerated beyond 
the threshold energy for spallation reactions to synthesize 
light elements. In our calculations, realistic massive star 
models are used as the initial conditions of SNe and the 
EOS takes into account the thermal radiation and ideal 
gas. We have compared the resultant energy distributions 
of ejecta with the empirical formula derived in Fields et 
al. (2002) for some SN explosions including furious and 
normal SNe and found that the energy distributions of 
ejecta from the numerical calculations and the empirical 
formula agree only in the high energy tail when the ex- 
plosion energy per unit ejecta mass significantly exceeds 
1.3 x 10 51 ergs/ M . Otherwise, the empirical formula 
overestimates or underestimates the ejecta mass at around 
the threshold energy for spallation reactions. Therefore it 
is necessary to numerically calculate SN explosions to ob- 
tain a correct energy distribution of ejecta for estimations 
of the yield of light elements. 
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Table 2 
Yields of light elements. 





b Li 


' Li 


y Bo 


iu B 


11 B 


log M Bc /M Q 


Model 


(10- 7 M ) 


(10- 7 M ) 


(10-7 M ) 


(10-7 M ) 


(lO- 7 M ) 




SN 1998bw 


2.38 


3.31 


0.999 


4.38 


13.4 


-8.0 


SN 2002ap 


0.0841 


0.114 


0.0348 


0.152 


0.464 


-8.8 


SN 19941 


0.0140 


0.0190 


0.00578 


0.0253 


0.0776 


-9.1 



To obtain the yields of light elements from the calcu- 
lated SN ejecta, we have numerically solved the transfer 
equation taking into account the energy loss due to ion- 
ization of the ISM and spallation reactions with the ISM. 
The results suggest that light elements synthesized from 
energetic SNe Ic like SN 1998bw by this mechanism can 
explain the enrichment of Be and B observed in the Galac- 
tic halo stars if the fraction of SN 1998bw like SNe in the 
early Galaxy was a factor of > 100 higher than current 
observational data suggest. This mechanism must not be 
the only primary mechanism that worked in the Galactic 
halo. Other primary mechanisms like superbubbles that 
supply light elements regardless of SN type are required to 
reproduce the observed abundance ratios such as Be/Fe. 
These SNe Ic can produce Li with more than one order 
of magnitude smaller amounts than indicated by observa- 
tions. However, lack of information on the Li abundances 
in the mctallicity range of — 2 £ [Fe/H] ,$ — 1 prevents us 
from deducing a firm conclusion. 

SNe are suggested to be associated with aspherical ex- 



plosions. The deviation from spherical symmetry will be 
able to increase the mass of ejecta with enough energies for 
spallation reactions for a given _B CX /M j because M(> e) 
is proportional to the ~ 3.4 — 3.6 power of this value. To 
illustrate this effect, a simplified situation will be consid- 
ered. Suppose the energy injected in the direction with 
a solid angle of n steradian is enhanced by a factor of 
two and the energy in the other directions is reduced by 
a factor of 1.5, then the empirical formula for the mass 
M(> e) indicates that this mass will increase by a factor 
of - 2 3 - 4 x 1/4 + (2/3) 3 - 4 x 3/4 - 2.8 while the total en- 
ergy will be unchanged. Applying the empirical formula 
obtained from the spherically symmetric calculations to 
this situation might lead to an erroneous result. Thus to 
further explore SNe Ic as a production site for light ele- 
ments, we need to perform multi-dimensional rclativistic 
hydrodynamic calculations for SNe Ic that can trace the 
motion of the outermost ejecta with a sufficient accuracy 
such as the calculations presented here. 
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